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ABSTRACT 



Within the Minimal Supersymmetric Grand Unified Theory (MSGUT) masses 
of the predicted supersymmetric particles are constrained by the world aver- 
aged values of the electroweak and strong coupling constants, the lower limits 
on the proton lifetime, the lower limit on the lifetime of the universe, which 
implies an upper limit on the dark matter density, the electroweak symmetry 
breaking originating from radiative corrections due to the heavy top quark, 
and the ratio of the masses of the b-quark and r-lepton. A combined fit shows 
that indeed the MSGUT model can satisfy all these constraints simultaneously 
and the corresponding values for all SUSY masses are given within the min- 
imal model, taking into account the complete second order renormalisation 
group equations for the couplings and the one-loop corrections to the Higgs 
potential for the calculation of the Mz mass and the Higgs masses. These one- 
loop corrections to Mz have been derived explicitly as function of the stop- 
and top masses and found to be small for the best solution, but unnaturally 
large for the 90% C.L. upper limits on the SUSY masses. 
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1. Introduction 

Grand Unified Theories (GUT) hold the promise of " explaining" the difference 
between the electromagnetic, weak and strong nuclear forces: their different strenghts 
are simply due to radiative corrections. Furthermore, they are candidates to explain 
several unrelated observations about our universe, e.g. they almost automatically lead 
to baryon number violation, thus providing a possible explanation for the matter- 
antimatter asymmetry in our universe [1] and the spontaneous symmetry breaking 
of the unified force into the known forces at a sufficient high energy can cause the 
inflationary scenario of the universe, thus providing an explanation for the origin of 
matter and the homogeneity of the universe on a large scale. The interested reader is 
refered to recent text books on this exciting subject [2] for more details and original 
references. 

The SU (5) group, which is the smallest group encompassing the SU (3) and 
SU(2) U(l) groups of the strong- and electroweak interactions, can be ruled out as 
a viable GUT, since it predicts too rapid proton decay. An upper limit on the proton 
lifetime can be estimated from the unification scale and any GUT is required to have 
the unification scale above fO 15 GeV in order to be compatible with the proton life 
time limits. This is not the case in the SU(5) model and these limits severely constrain 
other GUT models. 

Indeed, the coupling constants, as measured precisely at LEP, do not unify 
either within the SU (5) model, i.e. they do not become equal at a single energy if 
extrapolated to high energies, but within the supersymmetric extension of the SU (5) 
model (MSSM)[3], unification is obtained[4][5][6][7][8]. Supersymmetry[9] presup- 
poses a symmetry between fermions and bosons, thus introducing spin partners of 
the quarks and leptons - called squarks and sleptons - and spin 1/2 partners of the 
gauge bosons and Higgs particles - called gauginos and Higgsinos. Since these pre- 
dicted particles have not been observed sofar, these supersymmetric (SUSY) particles 
must be heavier than the known particles, implying that supersymmetry must be bro- 
ken. However, from the unification condition a first estimate of the SUSY breaking 
scale could be made: it was found to be of the order of 1000 GeV, or more precisely 
10 3±1 GeV [5]. The uncertainty in this scale is mainly caused by the uncertainty in 
the strong coupling constant. 

Clearly, the whole SUSY mass spectrum cannot be described by a single pa- 
rameter. In the MSSM one needs at least 5 parameters. So many parameters cannot 
be derived from the unification condition alone. However, further constraints can be 
considered: 

• M z predicted from electroweak symmetry breaking [10] [11]. 

• b-quark mass predicted from the unification of Yukawa couplings[12] [11] [13]. 

• Constraints from the lower limit on the proton lifetime[14] [15] [16]. 

• Constraints from the lower limit on the lifetime of the universe[17]. 

• The maximum value of the top mass is restricted by the couplings [11]. 
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• Experimental lower limits on SUSY masses[18]. 

All these constraints have been considered separately or partly [4] [5] [6] [7] [8] [15] [19] [20] . 
However, considering only one constraint at a time allows one to obtain only one re- 
lation between parameters. Trying to find complete solutions requires additional 
assumptions, like naturalness, no-scale models, fixed ratios for gaugino- and scalar 
masses or a fixed ratio for the Higgs mixing parameter and the scalar mass, assump- 
tions from supergravity, or combinations of these assumptions. 

The different assumptions lead sometimes to apparently conflicting results. 
For example, the dark matter constraint requires the mass of the scalar particles at 
Mqut to be below 500 GeV[17], while the analysis of the proton life time constraint 
requires this mass to be above 600 GeV[14]. 

It is the purpose of this paper to study all constraints simultaneously without 
any of the assumptions mentioned above in order to see if a solution within the 
minimal super symmetric model (MSSM) exists at all, which is non-trivial as should 
be clear from the contradictions mentioned before. 

It turns out that a solution exists indeed and the corresponding constraints 
on the SUSY mass spectrum, the strong coupling constant and the top mass are 
given, if all constraints from the experimental data mentioned above are imposed 
simultaneously (for the first time as fas as we know). 



2. Experimental constraints 

2.1. Unification of the couplings 

In the unified SU(2)(g)U(l) theory, the following well known tree- level rela- 
tions hold between the couplings and the gauge boson masses 

e = \fA-na = g sin 9w = g' cos 9w 
M w = \vg (2.1) 

M z = lvVg JT T^ 

from which it follows that 

p 2 a' 2 M 2 

^ 2 °w = - 2 = ^— 2 = l-^ (2-2) 
g 2 g' 2 + g 2 M 2 

Here g and g' are the couplings of the groups SU(2) and U(l) respectively, a is the 
fine structure constant and v is the vacuum expectation value of the Higgs field. 
If the model contains Higgs representations other than doublets, the theory has an 
additional degree of freedom, usually parametrized by the p parameter. 

In the SM based on the group SU(3) x SU(2) x U(l) we use the usual definitions 
of the couplings 

ai = (5/3)# /2 /(47r) = 5a/(3cos 2 9 w ) 

a 2 = # 2 /(47r) = a/sin 2 9 w (2.3) 

«3 = 9 2 slW) 
where g s is the SU(3) coupling. The factor of 5/3 in the definition of ot\ has been in- 
cluded for the proper normalization at the unification point[21]. The couplings, when 



4 



defined as effective values including loop corrections in the gauge boson propagators, 
become energy dependent ("running"). A running coupling requires the specification 
of a renormalization prescription, for which one usually uses the modified minimal 
subtraction (MS) scheme [22]. 

In this scheme the world averaged values of the couplings at the Z° energy are 

or\M z ) = 127.9 ±0.2 

sin 2 ^ = 0.2326 ±0.0006 (2.4) 
a 3 = 0.123 ±0.008 

The value of a' 1 is given in Ref. [23] and the value of sin 2 6-^$ has been presented 
at Marseille [24]. The 0:3 value corresponds to value measured at LEP form the ratio 
of the hadronic- and leptonic cross sections [25], which agrees well with the world 
average [26]. This value has the smallest theoretical uncertainty from the higher order 
corrections, we prefer to take this value alone. 

For SUSY models, the dimensional reduction DR scheme is a more appropriate 
renormalization scheme [27]. This scheme also has the advantage that all thresholds 
can be treated by simple step approximations. Thus unification occurs in the DR 
scheme if all three aj 1 ^) meet exactly at one point. This crossing point then gives 
the mass of the heavy gauge bosons. The MS and DR couplings differ by a small 
offset 

-J= = J— ^ (2.5) 
a DR a MS Un v 1 

where the C, are the quadratic Casimir coefficients of the group (C, = N for SU(iV) 
and for U(l) so ot\ stays the same). Throughout the following, we use the DR 
scheme for the MSSM. 

The energy dependence of the couplings is completely determined by the par- 
ticle content and their couplings inside the loop diagrams of the gauge bosons as 
expressed by the renormalization group (RG) equations. The RG equations can be 
rewritten as 

Af^ )= ^ + %& m+0{c ®)- (2 ' 6) 

In first order, i.e. all by = 0, the equations for the three a" 1 are independent with a 
linear solution in the a~ l — log /i plane. When the second order contributions are taken 
into account, the equations become coupled and the running of each a~ l depends on 
the values of the other two couplings. However, the second order effects are small 
because of the additional factor <x,/(47r) < 0.01. Higher orders are presumably even 
smaller by additional powers of <x,/(47r). We solve (2.6) by numerical integration. 

2.2. Mz constraint from electro weak symmetry breaking 

In the MSSM at least two Higgs doublets have to be introduced: 
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At the tree level the interactions of the Higgs fields can be parametrised by an 
effective potential of the form: 



2 fi 2 

V(H U H 2 ) = mlm 2 + m 2 2 \H 2 \ 2 - m\{E x E 2 + h.c.) + ^-[-^-(|#i| 2 - \H 2 \ 2 ) 2 , 

o 

where mi, m 2 and are mass parameters, for which we assume at the GUT 
scale the following boundary conditions: m\ = m\ = ji^ + ml, m\ = Bfi m . 

Radiative corrections from the heavy top and stop quarks can drive one of 
the Higgs masses negative, thus causing spontaneous symmetry breaking in the elec- 
troweak sector. In this case the Higgs potential does not have its minimum for all 
fields equal zero, but the the minimum is obtained for non-zero vacuum expectation 
values of the fields: 

v v 
< Hi >= vi = —j= cos (3, < H 2 >= v 2 = —j= sin (3, 
v2 v 2 

J = V 1±A^ tan^ 
2 v 2 

The scale, where symmetry breaking occurs depends on the starting values 
of the mass parameters at the GUT scale, the top mass and the evolution of the 
couplings and masses. This gives strong constraints between the known Z° mass and 
the SUSY mass parameters, as demonstrated e.g. in Ref.[10]. 

Minimization of the tree level potential yields: 

2 &{m\ - ml tan 2 (3) . 2m\ 



v= ,; / 9w , ' \ , sin 2(3 



(g 2 + g' 2 )(tan 2 (3 - 1)' m 2 + m 2 2 



M 2 = 9 -v 2 , M 2 = 9 -^-v 2 , 



4 ' * 4 

After including the one-loop corrections to the potential [28], the Z° mass becomes 
dependent on the top- and stop quark masses too. In this case we derived the following 
expression: 



M 2 2 m l -m 2 2 tan 2 p- A z 



Z to„2 



tan 2 (3 - 1 



A 



2 



3S 2 



,2 



z ~ 32tt 2 M 2 



f(m 2 tl ) + f(mj 2 ) + 2mf + (Aj - ^ 2 cot 2 (3) 



W L "Hi "H2 



f( m 2 tl )-f(m 2 t2 ) 
mfi — m 



where 



2 2 2 

el 1\ In ^ i\ i I 2\ ^ i ^ 

/ [m ) = m (In — - 1), h(m ) = — — In — 

mf m A — m z q m z q 

d(mlm 2 2 ) = 2- m l + m l \n^ 
m\ — m 2 m 2 
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rh 2 is the mass of a light squark and fh\ 2 the stop quark masses. Note that the 
corrections A z are zero if the top- and stop quark masses are identical, i.e. if su- 
persymmetry would be exact. They grow with the difference m 2 — m t 2 , so these 
corrections become unnaturally large for large values of the stop masses, as will be 
discussed later. 

The masses of the physical Higgs particles after spontaneous symmetry break- 
ing become after inclusion of the one-loop corrections [28]: 



m , 



A 2 
m 2 H± 
A 2 



A 2 
sin 2/3' 



3g 2 m 2 A t fi j \m 2 tl ) - j \m 2 t2 ) 



2ms IGttW^M 2 , 



= m 2 A + M^ + A 2 H , 



3^ 2 



m 2 /i 2 



h{m 2 n ) 



m 2 tl - m 2 t2 



h{m 2 2 ) 



32tt 2 sin 4 pMl 



w 



™h - m 2 t2 



m 



h,H 



m 2 A + M 2 z + An + A 



22 



±22 



^12 



i 



M 2 



An + A 22 ) 2 
4(cos 2 f3M 2 z + sin 2 f3M 2 A ) A 22 
4sin2/3(M 2 + M\)A 12 



-Am 2 A M 2 z cos 2 2(3 - 4(A n A 22 
-4(sin 2 (3M Z + cos 2 PM 2 A )A U 



A? 2 ) 



m. 



16tt 2 sin 2 f3M^ 



fx(A t + ficotfi) 



m t 2 i 



rn 2 2 



d(rh tl ,m t2 ), 



16tt 2 sin 2 f3Ml 



W L 



ln( 



m 2 tl m 2 t2 . 



+ 



2 A t {A t + n cot p) 
m 2 i - m 2 t2 



~2 



ln( 



in 



a ■ 



™t2 



A t (A t + n cot 0) 



Zr.2 



m 



™2 



tl 



m 



t2 



d(m tl ,rh t2 ) 



3g 2 



m\ /j,(A t + /j, cot P) 



16tt 2 sin 2 (3M 2 



w 



m 2 i 



m 2 2 



m 2 tl , A t (A t + n cot p) 



m 2 t2 



m 2 n 



™t2 



d{m 2 n ,rh 2 t2 ) 



A{t) = 



An 



1 + 6Y F(t) 



+ 



mi/2 
m 



Ho 



QY H 3 



l + 6YoF(t) t 



Here H 2 , H 3 , E and F are functions of the couplings as defined in [11]. 

2.3. Evolution of the masses 

In the soft breaking term of the Lagrangian m and rai/ 2 are the universal 
masses of the gauginos and scalar particles at the GUT scale, respectively and \i 
determines the masses of the particles in the Higgs sector. At lower energies the 
masses of the SUSY particles start to differ from these universal masses due to the 
radiative corrections. E.g. the coloured particles get contributions proportional to 
a 2 from gluon loops, while the non-coloured ones get contributions depending on 
the electroweak coupling constants only. The evolution of the masses is given by 



7 



the renormalisation group equations [11]. We have used analytical solutions[29] of 
these equations including the top Yukawa coupling, the Higgs mass parameters, mass 
mixing between the top quarks, mixing between neutralinos, mixing between the 
charginos and one loop radiative corrections to the Higgs potential. 

2.4- b- quark mass constraint 

Unification of the Yukawa couplings for a given generation at the GUT scale 
predicts relations for quark and lepton masses within a given family. This does not 
work for the light quarks, but the ratio of b-quark and r-lepton masses can be correctly 
predicted by the radiative corrections to the masses [11]. 

2.5. Proton lifetime constraints 

GUT's predict proton decay and the present lower limits on the proton lifetime 
yield quite strong constraints on the GUT scale and the SUSY parameters. As men- 
tioned at the beginning, the direct decay p — > e + 7r° via s-channel exchange requires 
the GUT scale to be above 10 15 GeV. This is not fulfilled in the SM, but always ful- 
filled in the MSSM. Therefore we do not consider this constraint. However, the decay 
via box diagrams with winos and Higgsinos predict much shorter lifetimes, especially 
in the preferred mode p — > VK + . From the present experimental lower limit of 10 32 yr 
for this decay mode Arnowitt and Nath[14] deduce an upper limit on the parameter 
B: 

B < 293±A2(M H J3M GUT ) GeV' 1 

Here M H3 is the Higgsino mass, which is expected to be of the order of M GUT . 
To obtain a conservative upper limit on B, we allow M H . A to become an order of 
magnitude heavier than Mqut, so we require 

B < 977 ± 140 GeV' 1 . 

The uncertainties from the unknown heavy Higgs mass are large compared with the 
contributions from the first and third generation, which contribute through the mix- 
ing in the CKM matrix. Therefore we only consider the second order generation 
contribution, which can be written as[14] : 

B = -2(a 2 /(a 3 sin(2/3)) (m 5 /m?) 10 6 

One observes that the upper limit on B favours small gluino masses rrig, large squark 
masses rriq, and small values of tan [3. To fulfill this constraint requires 

tan/3 < 5 

for the whole parameter space. Arnowitt and Nath note that requiring the gluino 
mass to be below 500 GeV implies the mass of the scalar particles (mo) at the GUT 
scale to be above 600 GeV. We will not impose this requirement on the gluino mass. 
Furthermore, they require M H , S < 3 M GUT , so they obtained tighter limits on tan (3, 
since we allow M H:i < 10 M GUT . 
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2.6. t-quark mass constraints 



For large Yukawa couplings the masses become dependent only on the cou- 
plings (pole term in the RGE). This yields a maximum value of the top mass, which 
has to be higher than the experimental mass, thus constraining the SUSY parameters. 
More precisely, the top mass can be expressed as: 

m t 2 = (An) 2 Y t (t) v 2 sm 2 ((3) 

M 2 

where the running of the Yukawa coupling as function of t — log(-^-) is given by[ll]: 

Y E(t) 



Yt(t) = 



1 + 6Y F(t) 



One observes that for large Yukawa couplings Y at the GUT scale, Y t becomes inde- 
pendent of Y and the maximum value of m t becomes: 



2 (4vr) 2 E{t) 



where E and F are functions of the couplings only[ll]. Clearly the experimental 
values of m t have to be below this upper bound, which is most easily fulfilled for 
larger values of tan/3. However, the proton life time limits require tan/3 < 5, as 
discussed before. In this case the upper limit on m t implies a constraint on the ratio 
of E/F, i.e. on the starting point at the GUT scale and the intermediate SUSY 
thresholds. 



2. 7. Constraints from the lifetime of the universe 

The lightest supersymmetric particle (LSP) is supposedly stable and would be 
an ideal candidate for dark matter. So the MSSM predicts dark matter. However, 
from the long lifetime of the universe one knows that the density of the universe cannot 
be higher than the critical density, which implies most of the LSP's have annihilated 
into photons. This can only happen fast enough if the squarks and sleptons are suffi- 
ciently light, thus posing a strong upper limit on some of the SUSY mass parameters. 
Requiring the density of the universe to be below the critical density translates into 
an upper bound of about 500 GeV on m for a large range of mi/ 2 [17], i.e. 

m < 500 GeV. 

2.8. Experimetal lower limits on SUSY masses 

SUSY particles have not been found sofar and from the searches at LEP one 
knows that the lower limit on the charged particles is about half the Z° mass (45 
GeV) and the Higgs mass has to be above 60 GeV[18]. This requires also minimal 
values for the SUSY mass parameters. 
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3. Constrained Fits 

3.1. Fit strategy 

As mentioned before, given the 5 parameters in the MSSM and olqut and 
Mqut, all other SUSY masses, the b-quark mass, M z , and the proton lifetime can 
be calculated. In addition, the complete evolution of the couplings including all 
thresholds can be performed. Furthermore, the dark matter constraint requires itlq 
to be below 500 GeV. 

Therefore we have adopted the following strategy: we varied m between 
and 500 GeV and fitted the remaining 5 parameters: aaur, Mqut, m>i/2, A*, and 
tan/3. The trilinear coupling Aq in the Higgs potential at Mqvt was kept mostly at 
0, but the large radiative corrections to it were taken into account for the top quarks. 
Varying A between +m and — m did not change the results significantly, so it was 
kept zero for the results quoted hereafter. 

The remaining parameters were fitted with MINUIT by minimizing the fol- 
lowing x 2 function: 



2 ^ {al\M z ) - g-^ SSMi {M z )) 2 

X 2^ 2 

1=1 °t 
(M z -91.18) 2 



(m b - 4.25) J 



(B -f' 7) W> 997) 



(D(mlm2m3)f 



(for D > 0) 



The first term is the contribution of the difference between the calculated and mea- 
sured coupling constants at Mz and the following three terms the contributions from 
the M^-mass, m^-mass, and m t -mass constraints. The last three terms give the con- 
tributions from the proton lifetime, the requirement of electroweak symmetry break- 
ing, i.e. D = m\m\ — m\ < 0, and experimental lower limits on the SUSY masses. 
The following errors were attributed: <7j equal the experimental errors in the coupling 
constants, as defined before, <r&=0.3 GeV, <tb=0.14 GeV, and all the other errors were 
set to 10 GeV. The value of these errors turned out not to be critical at all, since the 
corresponding terms in the numerator were usually zero in case of a good fit and even 
for the 90% C.L. values these constraints could be fulfilled and the x 2 was determined 
by the other terms, for which we know the errors. 
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For unification in the DR scheme, all three couplings a^ l (n) must cross at a 
single unification point in the a~ l — \i plane given by Mqut and «5ut (^ ne inverse 
of the unified coupling). Thus in these models we can fit the couplings at Mz by 
extrapolation from a single point at Mqut back to Mz for each of the aj's and taking 
into account all thresholds. Between the highest SUSY threshold and M z only the 
first order coefficients in the RGE are known for the individual thresholds, at least 
as far as we know. However, the second order coefficients must be between the 
values of the MSSM including all particles and the SM. So we varied these second 
order coefficients in this range for the small region between Mz and the highest SUSY 
mass. The difference was found to be negligible. Of course, for the large extrapolation 
between the highest SUSY mass and M GUT the complete second order RGE was used. 

3.2. Results 

The upper part of Fig. 1 shows the evolution of the coupling constants in the 
MSSM for two cases: one for the minimum value of the \ 2 function as defined above 
(solid lines) and one corresponding to the 90% C.L. upper limit of the thresholds of 
the light SUSY particles (dashed lines). The light thresholds are indicated in the 
lower function as the change in the first order coefficient in the f3 function, which 
corresponds to the first order change in the slopes of the curves at the top. One 
observes that the change in a 3 occurs in a rather narrow energy regime, corresponding 
to the threshold of squarks and gluinos, while for the other coupling constants the 
sleptons, higgsinos and winos contribute in addition, thus causing a somewhat more 
smeared threshold region. 

The parameters corresponding to these fits are tabulated in Table 1. The 
initial choices of m as well as the fitted paramters are shown at the bottom. As 
mentioned before, varying A between +m and — m does not influence the results 
very much, so it was kept at 0, but the large radiative corrections to it were taken into 
account. The remaining parameters do not depend strongly on the initial choices of 
mo and m t , as shown in Fig. 2. The value of mo was varied between and 500 GeV 
and m t between 140 and 170 GeV. Larger values of m are not allowed by the dark 
matter constraint and the range of m t is the preferred range from the b-quark mass, 
as shown in Fig. 3. The difference between the two bent lines originates from the 
difference in cigut and Mqut, as found from the fits for different input values of mj. 
The horizontal band corresponds to the mass of the b-quark after QCD corrections: 
nib = 4.25 ± 0.1 GeV[30]. Although the error of 0.1 GeV is the quoted one, we used 
the more conservative error estimate of 0.3 GeV, since the value of the mass depends 
on the value of the strong coupling constant at this mass, which is not too well known. 
Note that the allowed range of m t is in excellent agreement with LEP results too[18]. 

In order to obtain 90% C.L. upper limits on the parameters, we varied them 
until the x 2 value increased to 1.64. Of course, this depends on the errors, so we used 
the rather conservative errors defined above. The upper limit on tan (3 was found to 
be about 5; here the main constraint is coming from the proton lifetime. 

Since mo is rather small, the solutions are not strongly dependent on it, but 
they are mainly determined by \x and m^. Since also tan /3 is constrained to a rather 
narrow range, we are left effectively with two parameters, mi/ 2 and /i. However, they 
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were found to be strongly correlated, as shown if Fig. 4. In the minimum the x 2 value 
is zero, but one notices a long valley, where the x 2 is only slowly increasing. One 
can easily understand such a behaviour from Fig. 1: as long as the breakpoints in all 
curves move up simultaneously, unification can be obtained at a higher value of olgut, 
as is obvious from a comparison of the dashed and solid lines. If one tries to combine 
the solid lines with the dashed lines into a single unification point, this clearly does 
not work. Consequently the breakpoints in all curves have to be close together in 
energy. However, a 3 is independent of \x in contrast to the breakpoints in ct\ and ct 2 - 
This implies that both ji and have to increase simultaneously in order to keep 
the breakpoints in all three curves close together. This strong correlation has been 
neglected in previous analysis, where ji and were choosen independently [10] [31]. 
This clearly cannot work in a more detailed analysis, as observed too in Ref. [32]. 

The steep walls in Fig. 4 originate from the experimental lower limits on the 
SUSY masses and the value of Mz from radiative symmetry breaking. 

Fig. 5 shows the evolution of the masses for the minimum value of x 2 an d 
m =500 GeV. From Table 1 one observes that some SUSY masses can go up to several 
TeV, if one considers the 90% C.L. Such large values spoil the cancellation of the 
quadratic divergencies. If supersymmetry would be exact, i.e. as long as the masses 
and couplings of the particles and their superpartners are the same, the contributions 
of fermions and bosons in the loops would exactly cancel each other, thus eliminating 
in first order all divergencies. This can be seen explicitly in the corrections to M z : 
Az is exactly zero if the masses of stop- and top quarks are identical. For the SUSY 
masses at the minimum value of x 2 the corrections to Mz are indeed very small, as 
shown in Fig. 6 on the left hand side, but for the solutions corresponding to the 90% 
C.L. upper limit the corrections to M z are about 50 times M z itself (see Fig. 6 right 
hand side). In addition, the value of the determinant l l A ={m\m 2 2 —ml) 1 / 4 in the Higgs 
potential times its sign is shown as a dotted line. It clearly becomes negative, which 
implies the potential takes the shape of the mexican hat. The evolution of the masses 
in the Higgs potential is shown in the lower parts of Fig. 6. 

If we require that only solutions are allowed, for which the corrections to M z 
are not large compared with M z itself, we have to limit the mass of the heaviest stop 
quark to about 1 TeV. In this case the 90% C.L. upper limits if the individual SUSY 
particles are given in the right hand column if Table 1. The correction to Mz = 5 
times M z in this case. 

To determine the lower limits is somewhat more cumbersome, since not all 
particles reach their minimum for the same set of parameters. We used the following 
strategy: each parameter was scanned for the lowest possible value, while keeping the 
other parameters free. A zero value for the scalar mass m could not be excluded, 
which corresponds to the so-called no-scale model and the corresponding results have 
been given in the left column of Table 1. The parameters /i and m\/2 cannot be 
choosen independently and their minimal values were obtained for m = 400 GeV, as 
shown in the second column of Table 1. Lower values would cause a wino mass below 
the experimental lower limit of about 45 GeV. The top mass for these minima was 
kept at 170 GeV, since smaller values would give a too high value of m h . 

One observes from Table 1 the well known effect [2 8] that the Higgs particle, 
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which gets a negative mass squared before spontaneous symmetry breaking (SSB), 
gets only a small mass after SSB. The mass of this particle, called h in Table 1, is a 
rather strong function of m t , as shown in Fig. 7. For each value of m t , the parameters 
TO1/2, iii and tan/5 were determined by minimising the x 2 f° r the indicated value of 
m . One observes that the mass of the lightest Higgs particle varies between 60 and 
120 GeV. These values correspond to the minimal value of x 2 , but even if the 90% 
C.L. limits are taken, the mass increases only to 125 GeV. 

4. Summary 

The MSSM model has many predictions, which can be compared with experi- 
ment, even in the energy range where the predicted SUSY particles are out of reach. 
Among these predictions: 

• M z . 

• m b . 

• Proton decay. 

• Dark Matter. 

• Upper limit on m t . 

It is surprising, that the minimal supersymmetric model can fulfil all experimen- 
tal constraints for these predictions. As far as we know, supersymmetric models 
are the only ones, which are consistent with all these observations simultaneously. 
Other models can yield unification too [6], but they do not exhibit the elegant sym- 
metry properties of supersymmetry, they offer no explanation for dark matter and 
no explanation for the electroweak symmetry breaking. Furthermore the quadratic 
divergencies do not cancel. 

From the above constraints we find at the 90% C.L. (see Table 1): 






< 


m 


< 


500 GeV 


70 


< 


mi/2 


< 


3410 (475) GeV 


500 


< 


n 


< 


3150 (1000) GeV 


1 


< 


tan/? 


< 


5 


138 


< 


m t 


< 


186 GeV 


0.108 


< 


a s 


< 


0.132 



The upper limit on m originates from the dark matter constraint, the upper 
limit on tan/3 from the lower limit on the proton lifetime. The fact that tan/5 is 
so much smaller than the ratio of top- and b-quark mass implies that the Yukawa 
coupling of the b-quark is negligibly small, so one does not have to consider its 
contributions in the renormalisation group equations. 



13 



Good fits are only obtained for a s between 0.108 and 0.132, if the error on a s 
is taken to be 0.008. The bottom mass constraint together with the given couplings 
require the top mass to be between 138 and 186 GeV. 

The values in brackets indicate the 90% C.L., if one requires that the exact 
one-loop corrections to Mz are not large compared with Mz itself, which requires the 
heaviest stop quark to be below 1000 GeV. In this case the correction to M z is about 
5 M^and the corresponding constraints on the other SUSY masses are (see Table 1 
for details): 



25 


< 


7(X?) 


< 


200 GeV 


45 


< 




< 


385 GeV 


190 


< 


9 


< 


1110 GeV 


430 


< 


Q 


< 


1075 GeV 


200 


< 


ti 


< 


730 GeV 


390 


< 




< 


1000 GeV 


235 


< 


e L 


< 


520 GeV 


130 


< 


£r 


< 


440 GeV 


335 


< 


H 


< 


785 GeV 


525 


< 


H 


< 


870 GeV 


60 


< 


h 


< 


125 GeV 



The values of m 1 / 2 and ji are positively correlated, i.e. a large (small) value 
of mi/2 corresponds to a large (small) value of fx, as is apparent from Fig. 4. This 
strong correlation was usually not taken into account in previous analysis, in which 
mi/2 and /i were restricted by ad-hoc assumptions [10] [31]. 

The detailed mass spectra have been given in Table 1. The lightest Higgs 
particle is certainly within reach of experiments at present or future accelerators. Its 
observation in the predicted mass range of 60 to 125 GeV would be a strong case in 
support of this minimal version of the supersymmetric grand unified theory. 
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particles 


masses in [GeV] 


Lower limit 


best fit 


Upper limit 90% C.L. 


fine tuning 


no 


no 


no 


no 


mi < 1 TeV 


/ VAl/ 


134 


23 


46 


1501 


199 


IA2/ 


254 


46 


88 


2850 


383 


W(yf) 

VAl / 


254 


44 


87 


2900 


385 


9 


799 


193 


330 


6826 


1109 


e t 


236 


406 


412 


2331 


521 


e n 


131 


402 


404 


1384 


441 


VL 


228 


399 


406 


2331 


516 




725 


434 


497 


6035 


1075 


Or? 


697 


431 


490 


5700 


1036 




498 


201 


226 


4740 


729 




748 


391 


472 


5600 


1000 




550 


321 


417 


3311 


771 


AA4/ 


572 


337 


436 


3323 


784 




569 


337 


433 


3324 


783 




91 


101 


96 


115 


127 


H 


624 


525 


636 


2897 


871 


A 


619 


523 


633 


2897 


870 


H± 


625 


529 


638 


2897 


873 


SUSY parameters 


m 





400 


400 


500 


400 


mi/2 


329 


70 


172 


3413 


475 


n 


550 


507 


576 


3150 


1009 


tan j3 


2.0 


3.7 


2.2 


2.9 


3.5 


m t 


147 


186 


172 


138 


175 




24.7 


24.2 


24.3 


26.4 


25.1 


M GUT 


1.3 10 16 


2.0 10 16 


1.9 10 16 


0.7 10 16 


1.2 10 16 



Table 1: Values of SUSY masses and parameters. For the lower limits of the SUSY 
masses one should take the lowest value from the first two columns, except for the 
value of the lightest Higgs h. In that case the lowest value is given by Fig. 7. 
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Figure 1: Evolution of the inverse of the three couplings in the MSSM. The line 
above Mqut follows the prediction from the supersymmetric SU(5) model. The SUSY 
thresholds have been indicated in the lower part of the curve: they are treated as step 
functions in the first order f3 coefficients in the renormalisation group equations, which 
correspond to a change in slope in the evolution of the couplings in the top figure. 
The dashed lines correspond to the 90% C.L. upper limit for the SUSY thresholds. 



18 



> 

O 



> 

CD 

o 



CO. 



2500 
2000 
1500 
1000 
500 




10 
8 
6 
4 
2 




2 
1.5 
1 

0.5 
















K)G< 
'0G( 


















m top- 14 
m top- 1 


?v- 
















































1 1 1 


i i i i i i i i i i i i i i i i i i i i 1 i i i i 1 i i i i T i 1 ~i T 


i ' i i i 


' i i i i 1 i i 



50 100 150 200 250 300 350 400 450 500 

m [GeV] 




50 100 150 200 250 300 350 400 450 500 

m [GeV] 



50 100 150 200 250 300 350 400 450 500 

m [GeV] 



50 100 150 200 250 300 350 400 450 500 

m [GeV] 



Figure 2: The fitted MSSM parameters as function of the choice of m : above mi/ 2 , 
middle ji and below tan/3. The shaded area indicates the variation, if m t is varied 
between 140 and 170 GeV. 
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Figure 3: The correlation between m& and m t for mo = 500 GeV. The hatched 
area indicates the experimental value for rrib = 4.25 ± 0.3 GeV (the crossed area for 
mi, = 4.25 ± 0.1 GeV) and the curved lines give the correlation for different values of 
the coupling constant at the GUT scale. 
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Figure 4: The correlation between m\ji and \i for m =500 GeV. The vertical axis 
gives the value of x 2 - 
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Figure 5: The evolution of masses of the SUSY particles. The lines indicated as 
photino and zino are the lines corresponding to the two lightest neutralinos, while 
the lines of wino and charged higgsino correspond to the two lightest charginos. The 
lightest Higgs (neutral Higgs 2 in lower left corner) gets its small mass only after 
spontaneous symmetry breaking via the Higgs mechanism. 
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Figure 6: The evolution of the Mz mass, the determinant in the Higgs potential, 
and the value of the Higgs mass squared for the optimum solution (left) and the 
solution corresponding to the 90% C.L. upper limits (right). The dashed lines in 
the top curves show the solution of the minimization of the Higgs potential at the 
tree level. This line coincides with the solid line for the best solution, as expected if 
supersymmetry is not badly broken, but the one-loop corrections to the Mz masses 
become unnaturally large for the right hand side. 
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Figure 7: The mass of the lightest Higgs particle as function of the top quark mass. 
The parameters of tan/5, mo, and mi/ 2 were optimised for each choice of m t and m . 
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